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Abstract. - We present the results of electronic structure calculations for the two-dimensional 
ferromagnet Rb2CrCl4. They are obtained by the augmented spherical wave method as based 
on density functional theory and the local density approximation. In agreement with exper- 
imental data Rb2CrCl4 is found to be semiconducting and displays long-range ferromagnetic 
order of the localized Cr 3d moments. The magnetic properties are almost independent of the 
structural modifications arising from the Jahn- Teller instability, which leads from the parent 
body-centered tetragonal K2NiF4 structure to a side-centered orthorhombic lattice. In contrast, 
our calculations give evidence for a strong response of the optical band gap to the corresponding 
structural changes. 



Introduction. ~ Since the discovery of CrBrs as a ferromagnetic semiconductor by Tsub- 
okawa [1] this narrow class of materials has attracted increasing interest. To some part this 
is for fundamental reasons, since these compounds serve as ideal candidates for the Heisen- 
berg model. However, the prospect of using ferromagnetic semiconductors in spin electronics 
has motivated an even greater amount of research activities in the last years. Only recently, 
La2MnNi06 has been shown to exhibit long-range ferromagnetic order below Tc ~ 280 K 
with a magnetization of 5 /ib per formula unit [2] and an optical band gap of about 0.6 eV as 
determined by electronic structure calculations [3]. 

Rubidium tetrachlorochromide, Rb2CrCl4, also belongs to this class of non-metallic fer- 
romagnets albeit with a much lower ordering temperature of Tc = 52.4 K [4,5]. According 
to neutron scattering measurements of the spin-wave dispersion this compound behaves as 
a two-dimensional, S* = 2, easy-plane ferromagnet with an in-plane exchange constant of 
J = 7.7 K [6]. Since the interlayer coupling was estimated to be about four orders of mag- 
nitude smaller than the intralayer interaction, three-dimensional long-range order was also 
attributed to the weak uniaxial anisotropy. The latter, in addition, leads to a small energy 
gap of about 0.1 meV in the spin wave dispersion at the center of the Brillouin zone, which 
precludes true XY-type behaviour [5,7]. Furthermore, the anisotropy causes a canting of the 
ordered spins of 5.2° away from the (110) direction [7]. From susceptibility measurements 
and neutron scattering investigations of the static critical properties a crossover from an ideal 
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two-dimensional XY-likc behaviour to a narrow three-dimensional region on approaching Tc 
from either side of the transition was inferred [8-10]. 

The crystal structure of Rb2CrCl4 is closely related to the body-centered tetragonal struc- 
ture of K2NiF4 with space group D\J^ (I4/mmm). In this structure the magnetic ions are 
arranged on a square planar lattice and coupled via the anions located midway between them. 
However, the Cr electronic configuration in Rb2CrCl4 with triply occupied t2g states and an 
orbital degeneracy in the singly occupied Cg states is susceptible to a Jahn- Teller instability. 
Indeed, the situation is not unlike that foimd for K2CUF4, which has a Cu (P configuration. 
While for the copper compound early experiments were interpreted in favour of the K2NiF4 
structure with octahedra compressed parallel to the c-axis, Khomskii and Kugel argued on 
the basis of theoretical considerations that this structure was incompatible with the observed 
ferromagnetic order [11,12]. Subsequently performed neutron and x-ray diffraction studies 
confirmed these ideas and gave evidence for a side-centered orthorhombic lattice with space 
group D^fi (Cmca) and two formula units per cell [13,14]. In this structure, the octahedra 
are elongated with the long axis lying in the plane and pointing alternatively parallel to the 
a- and 6-axis of the parent K2NiF4 unit cell. For Rb2CrCl4, first evidence for the orthorhom- 
bic structure came from NMR experiments, which revealed that the local environments of 
the chromium sites are elongated chlorine octahedra having almost tetragonal symmetry [15]. 
Lateron, the space group D^f^ (Cmca) was also confirmed for Rb2CrCl4 [7, 16, 17]. 

In this space group the Rb atoms occupy the Wyckoff positions (8d), (0,0, ±zp{b), with 
ZRh = 0.3568 and the Cr atoms are located at the (4a) sites at (0, 0, 0). The apical CI atoms are 
found at the Wyckoff positions (8d) with zci = 0.1508, whereas the equatorial CI atoms occupy 
the (8f) positions at (1/4 + (5, 1/4 + 5, 0). Here, the distortion parameter S characterizes the 
deviation from the ideal K2NiF4-strueture; the experimental value of (5 = 0.0154 corresponds 
to a 0.16 A shift of the CI ions off the central positions between two neighbouring Cr sites [7] . 
As a consequence, the in-plane Cr-Cl bond distances amount to 2.70 and 2.39 A, whereas the 
apical distance is 2.37 A. 

Despite their analogous crystal structures there is one important difference distinguish- 
ing Rb2CrCl4 from K2CUF4. In the latter compound, the single unpaired electron in the 
configuration of the Cu^+ ion occupies the dz^_x2 or d^2_j,2 orbital, depending on the ori- 
entation of the local octahedron. Due to the peculiarities of the crystal structure the local 
axes of neighbouring Cu sites are orthogonal and so are the orbitals. Hence, according to the 
Goodenough-Kanamori-Anderson (GKA) rules ferromagnetic order results. In contrast, the 

configuration of the Cr^+ ions is characterized by a single Cg electron in a d^yi_^2 or d^^2_^2 
orbital. Yet, despite their 90° rotation, orbitals of this kind placed at neighbouring sites are 
not orthogonal and the GKA rules would be in favour of an antiferromagnetic coupling. 

In order to resolve this issue in a more systematic manner, Feldkemper and Weber started 
from a multiband Hubbard model, which included both the metal d and the ligand p or- 
bitals [18]. Using a Rayleigh-Schrodinger perturbation expansion to map this model onto 
an effective Heisenberg Hamiltonian they were able to describe the superexchange coupling 
between magnetic sites. Due to the inclusion of various interaction paths their calculations 
extended the basic GKA rules and resulted in a net ferromagnetic interplanar coupling in 
agreement with the experimental findings. 

In this Letter, we take account of the obvious lack of first principles studies for Rb2CrCl4 
and complement the model approaches by electronic structure calculations, thereby continuing 
previous work on ferromagnetic semiconductors [19-22]. In doing so, we aim especially at ex- 
tending the current understanding of the interplay between the Jahn- Teller distortion and the 
optical band gap as well as the magnetic ordering. Previous electronic structure calculations 
for K2CUF4, which confirmed the non-metallic ferromagnetic ground state, demonstrated. 
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that neither the lattice distortion due to the Jahn- Teller instability nor the magnetic order 
alone are sufficient to produce the optical band gap [20]. Since ferromagnetic ordering was 
also found for the ideal K2NiF4 structure {S ~ 0), it is the semiconducting ground state, which 
requires the Jahn- Teller effect, and not the magnetic order. Our results for Rb2CrCl4 are in 
line with this interpretation. 

Methodology. - The calculations are based on density functional theory and the local 
density approximation as implemented in the scalar-relativistic augmented spherical wave 
method [23,24]. In order to represent the correct shape of the crystal potential in the large 
voids of the open crystal structure, additional augmentation spheres were inserted. Optimal 
augmentation sphere positions as well as radii of all spheres were automatically generated by 
the sphere geometry optimization (SGO) algorithm [25]. Self-consistency was achieved by an 
efficient algorithm for convergence acceleration [26]. Brillouin zone sampling was done using 
an increasing number of k-points ranging up to 576 points within the irreducible wedge. 

Results and Discussion. - Partial densities of states as grown out of both spin-degenerate 
and spin-polarized calculations are displayed in Fig. ^ Of the three groups of bands the first 
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Fig. 1 - Partial Cr 3ci and CI 2p DOS for the experimental crystal structure [5 = 0.0154). Here and 
in the following figures results of spin-degenerate and spin-polarized calculations axe given on the left 
and right, respectively. 

one extending from about — 6.5eV to — 3eV originates predominantly from the CI 3p states. 
In contrast, the second and third group as found at and above the Fermi level trace back 
mainly to the Cr id t2g and eg states, d-p hybridization, leading to admixtures of these 
states in the energy region, where the respective other partner dominates, is much reduced. 
Finally, taking into account spin-polarization by allowing for long-range ferromagnetic order, 
we observe an almost rigid exchange splitting of the d states of « 2.8 eV going along with 
a total energy lowering of about 0.3 Ryd per unit cell. The calculated magnetic moment of 
3.718/iB per chromium site agrees perfectly well with the value of 3.7±0.4/iB as deduced from 
neutron diffraction data [27]. 

Concentrating on the small energy region of the Cr 3d states we display in Figs. |5] and 
131 the electronic band structure and the partial Cr 3d densities of states. For the former 
we use the Brillouin zone of the side-centered orthorhombic lattice with the high-symmetry 
points Z = (1/2,0,1/2), X = (1/2,0,0), P = (1/2,1/2,0), and N = (0,1/2,0) (see Figs. 4 
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Fig. 2 - Electronic structure for the experimental crystal structure {S = 0.0154). In the spin-polarized 
calculation spin-majority and spin-minority bands are given in green and blue, respectively. 



and 5 of Ref. [20] for more details on the relation between the Brillouin zones of K2NiF4 and 
K2CuF4/Rb2CrCl4). The almost vanishing dispersion along the line T-Z reflects the two- 
dimensional nature of the crystal structure. In the representation of the partial densities of 
states we use the local coordinate system with the z axis pointing parallel to the elongated 
distance of the local CrCle octahedron. We distinguish the narrow Cr 3d t2g states from the 
d3-,2_^2 and dx^-y'^ contributions. Due to the elongation of the octahedra the d3j,2_^2 states 
are shifted to lower energies as compared to the d.j.2_y2 states. As a consequence, we observe 
energetical overlap of the d^^'i-r^ states with the t2g bands and a gap of about 0.5 eV between 
the eg states. 

While the spin-degenerate calculations lead to metallic behaviour, the exchange splitting 
places the Fermi energy between the spin-majority d^^'^^^'^ and d.j.2_y2 states, leaving an 
optical band gap of 0.56 eV. Our results thus clearly confirm the orbital ordering as deduced 
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Fig. 3 - Partial Cr 3rf DOS for the experimental crystal structure {S = 0.0154). The t2g contribution 
is scaled to one orbital. 
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from the polarized neutron diffraction data by Miinninghoff et al. with all five d orbitals singly 
occupied except for the dx2_y2 orbital, which is empty [28]. 

In addition to the above calculations we have considered an antiferromagnetic order with 
antiparallel spin alignment of the two Cr sites inside the side-centered orthorhombic unit cell. 
Again, the bands experience an almost rigid spin splitting, which gives room for an optical gap 
of 0.43 eV. The magnetic moments per chromium site amount to ±3.627/iB. Yet, the energy 
of the antiferromagnetic solution is by 27meV higher than that found for the ferromagnetic 
situation. 

All calculations presented so far were based on the experimental crystal structure, i.e. for 
5 = 0.0154. In order to study the effect of the distortion on the electronic and magnetic proper- 
ties we performed additional calculations for several values of 5 varying from zero to 1.5 times 
the experimental value. The results are summarized in TableQ] which comprises spin-polarized 
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Table I - Total energies (in Ryd per 2 Cr), magnetic moments (in /is), and optical band gap (in eV) 
for different distortion parameters 5. The experimental crystal structure corresponds to 5 = 0.0154. 

calculations with both ferro- and antiferromagnetic order. Obviously, the total energies for 
both situations assume a minimum for the experimental value of 5. However, for the antifer- 
romagnetic solutions total energies are by w 27meV above those for the ferromagnetic case. 
In contrast, the optical band gap for both kinds of spin-polarized calculations assumes a max- 
imum value at (5 = 0.0154; it vanishes for ferromagnetic order in the ideal K2NiF4 structure 
as has been also observed in our previous work on K2CUF4 [20]. Surprisingly, the magnetic 
moments at the Cr sites are almost independent of the distortion parameter 5, i.e. of the Cr-Cl 
bond lengths. This is due to the localized nature of these moments, which forms the basis 
for the treatment within the Heisenberg model, and reflects the finding of Feldkemper and 
Weber that a reasonable variation of hopping parameters changes the magnetic coupling only 
slightly [18]. 

Finally, we aim at understanding the role of the Jahn- Teller effect for the changes of the 
electronic states. To this end, we display band structures calculated for hypothetical crystal 
structures with zero distortion parameter 5 in Fig. 0] On comparing these results to Fig. [21 
we recognize as the main effect of the Jahn- Teller distortion the splitting of the bands in the 
middle of the line P-F as well as near the X- and N-points. While in the spin-degenerate 
case this splitting occurs at about 1 eV above the Fermi energy, it opens the optical band gap 
for the ferromagnetically ordered state. Surprisingly, the spin-majority bands for (5 = and 
finite distortion as displayed in these figures look very similar to the spin-degenerate bands 
of K2NiF4, where the band degeneracies occur exactly at the X- and N-points and thus give 
rise to perfect Fermi-surface nesting [20]. While the nesting is not complete here, one would 
nevertheless be tempted to interpret the Jahn- Teller effect as resulting from a Fermi-surface 
instability of the spin-majority bands. However, we point out that, due to the comparatively 
low Curie temperature of Rb2CrCl4, we cannot use the spin-polarized calculation as a starting 
point but would rather need the band structure of a paramagnetic system, which is not 
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available. 

Conclusion. - In conclusion, wc have reported on the first density functional based elec- 
tronic structure calculation for the ferromagnetic semiconductor Rb2CrCl4. In agreement 
with the available experimental data we obtain a two-dimensional disperion of the electronic 
bands and a non-metallic ferromagnetically ordered ground state. The well localized magnetic 
moments of 3.7 fiB are carried by the Cr 3d <2g and (13^2 orbitals, in perfect agreement 
with the magnetic moments and orbital occupations deduced from neutron diffraction ex- 
periments. Studying the influence of the Jahn- Teller distortion we find only minor influence 
on the size and ferromagnetic order of the magnetic moments, thus confirming the model 
approach by Feldkemper and Weber. In contrast, the electronic band structure experiences 
drastic changes on switching on the Jahn- Teller distortion turning, in particular, from metallic 
to semiconducting. 

* * * 
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